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Abstract 

We investigate the properties of the Dirac operator on manifolds with boundaries 
in presence of the Atiyah-Patodi-Singer boundary condition. An exact counting of the 
number of edge states for boundaries with isometry of a sphere is given. We show that the 
problem with the above boundary condition can be mapped to one where the manifold 
is extended beyond the boundary and the boundary condition is replaced by a delta 
function potential of suitable strength. We also briefly highlight how the problem of the 
self-adjointness of the operators in the presence of moving boundaries can be simplified 
by suitable transformations which render the boundary fixed and modify the Hamiltonian 
and the boundary condition to reflect the effect of moving boundary. 


1 Introduction 

Wolfgang Pauli once remarked ‘God made the bulk, but the devil invented the surfaces’ p]. 
The truth of the statement cannot be over emphasized since the presence of boundaries and 
boundary conditions makes the analysis of a problem much more non-trivial. Manifolds with 
boundaries are ubiquitous in physics - from topological insulators [2j and quantum Hall sys¬ 
tems 13 U in condensed matter physics to black hole horizons and the cosmological horizon in 
general relativity I3i- The dynamics of the system of interest is usually encoded in the form 
of suitable differential operators on these manifolds. An analysis of the differential operators 
on manifolds with boundaries then becomes essential. One of the most important questions 
in this context is that of self-adjointness of these operators. The two most thoroughly studied 
operators in mathematical physics are the Laplacian and the Dirac operator [7j. 

While a general analysis of the question of self-adjointness of these operators on manifolds 
with boundaries is possible, see mm for instance, it seems that an investigation of these issues 
in the presence of certain isometries of the underlying manifold might be physically useful 
since, in certain cases, the presence of these isometries either allows for an exact solution or 
provides sufficient simplification so that most of the relevant aspects of the problem can be 
understood even without having an exact solution. This is especially so when one considers 
systems of relevance in physics where often one is interested in understanding the qualitative 
behavior of the system under consideration or where the underlying geometry of the problem 
is simple enough to work out the quantitative details. Following this reasoning, in this paper, 
we focus on the properties of the Dirac operator on manifolds with boundaries. 

It is well known that the question of self-adjointness of operators is intimately related 
to that of the domain on which these operators are defined which, in turn, is related to the 
choice of the boundary conditions. The two most commonly used boundary conditions are the 
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Dirichlet and the Neumann boundary conditions for which either the function or its normal 
derivative vanishes at the boundary respectively. However, for the Laplacian, a more general 
class of boundary condition, the so called Robin boundary condition, is also possible. For 
the Robin boundary condition, instead of requiring the function or its normal derivative to 
vanish, one has the condition that the function and its normal derivative are related at the 
boundary: ip\dM = Kip'ldM- 

It is clear that the Dirichlet and the Neumann boundary conditions are two extreme 
cases of the Robin boundary condition corresponding to k = 0 and k = oo, respectively. 
More interestingly, we find that the Robin boundary condition brings in a non-trivial length 
scale in the problem since the parameter k has the dimensions of length [L], Some of the 
consequences of having this non-trivial length parameter were explored in [9] where it was 
shown that Robin boundary condition leads to the presence of edge states on the boundary, 
with the number of these edge states being related to k and to the geometrical information 
of the boundary, k plays the role of the penetration depth. 

The boundary conditions for the Dirac operator were given by the Atiyah, Patodi and 
Singer [TOl [TTi, 12]. A further generalization of the Atiyah-Patodi-Singer (APS) boundary 
condition was recently given in m- As will be seen in more detail subsequently, like the 
Robin boundary condition for the Laplacian, the generalized APS boundary condition brings 
in a length parameter (jit) in the problem. In principle, more general boundary conditions are 
possible but in this paper we will not analyse this issue further since the associated isometries 
are limited. 

In this paper we focus on the APS boundary condition and its generalization with examples 
where the underlying manifold and its boundary have the isometries of sphere. We show that 
these boundary conditions can be thought of as a generalization of the Robin boundary 
condition in the context of the Dirac operator. We also show that this boundary condition 
can lead to several edge states on the boundary whose precise number depends on the mass 
parameter and the geometry of the boundary. We give an analytic expression for the 

number of edge states as a function of these parameters. 

In real systems the boundaries are imposed by applying suitable potentials. We would like 
to incorporate this aspect even at the mathematical level and as a second important aspect of 
this paper we show how the effect of the boundary conditions can be mimicked by replacing the 
boundary (and with it the corresponding boundary condition) by a delta function potential 
at the location of the boundary. In the context of the Schroedinger equation in one dimension 
it is well known that a negative delta function potential can lead to the existence of a bound 
state which is generally achieved by demanding the continuity of the wave function across the 
5-function potential (so that the derivative of the wave function is discontinuous). However, 
as has been shown in m , one can incorporate not only a 5 term but also a 5' term in the 
potential. Unlike the usual case, the analysis proceeds more naturally by considering the 
continuity equation at the location of the 5, 5' terms. Interestingly, this can be generalized to 
higher dimensions also. 

We will show that incorporation of the (^-potential is> i n fact, required when working 
with the Dirac operator. As an added advantage of replacing the boundary condition by 
interactions we note that since there are no boundaries in the problem any more, the question 
of self-adjointness of the differential operator becomes simpler. We show that the role of the 
parameter //, which in the absence of the boundary condition is not present in the problem, 
is taken up by the strength of the delta function potential. We give an exact counting of the 
number of edge states for this problem and show how this is related to the number of edge 
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states in the presence of the generalized APS boundary condition. 

Here we would like to mention some of the recent works incorporating a 5 (and/or a 
S' potential). The effect of singular potentials on vacuum energy has been discussed in 
m ns ini m [T9] (including the interaction between extended objects [20]). These pa¬ 
pers discuss interesting situations where the boundary conditions and self-adjointness play 
important roles which can possibly be tested experimentally. In addition, the effect of dis¬ 
continuous potentials on heat kernel expansion has also received some attention [SO 112 with 
the focus on self-adjoint Laplacian operators which are positive (since the primary concern 
was QFTs which are stable). This puts the constraint of looking for boundary conditions 
satisfying this requirement. Later we have pointed out that stability can be considered at 
higher temperature or with i/> 4 interactions (like spontaneous symmetry breakdown) and an 
analysis in these circumstances for heat kernel will complement their analysis. 

Finally, to illustrate the significance of the delta function potential further, we show 
that the effect of time dependent (or moving) boundaries and boundary conditions can also 
be mimicked by a suitable delta function potential. Time dependent boundaries which are 
ubiquitous in physics, lead to time dependent domains in the Hilbert space. In this case 
the meaning of time derivative of the state requires careful analysis. Replacing the boundary 
with 5 and 5' potential enables physical meaning of these considerations. The effect of moving 
boundaries has been investigated briefly in the literature, see E2 122 GEE Eg E7], and these 
are of relevance in diverse physical situations like an evolving black hole horizon or the time 
dependent cosmological horizon in an expanding universe, the Casimir effect in the presence 
of moving boundaries to name a few. 

In short, the objective of this paper is to show that - (i) the APS boundary condition can 
be seen as generalized Robin boundary condition and the same method as used to find the 
number of edge states for the Laplacian in the presence of the latter [9] can be used to find 
the number of edge states for the Dirac operator in presence of the APS boundary condition, 
(ii) the effect of the Robin or the APS boundary conditions (as far as the edge states are 
concerned) can be mimicked by removing the boundary and putting a delta potential of 
suitable strength at the location of the boundary and (iii) time varying boundaries can be 
analyzed by considering fixed boundaries but with a time dependent strength of the potential. 

The paper is organized as follows. In Sec. 2 after briefly reviewing the introduction 
of generalized APS boundary condition, we consider examples of euclidean M 2 and M 3 with 
a two-disc T> 2 and a three-ball £> :i removed. Imposing the the generalized APS boundary 
condition we find the exact number of edge states as a function of the parameters in the 
problem. Next, in Sec. 3, we show how the effect of the boundary condition can be mimicked 
by introducing a delta function potential of suitable strength at the location of the boundary. 
In Sec. 4 we further illustrate the utility of the delta function potential by showing how it can 
be used to study the effects of moving boundaries. In Sec. 5 we discuss the possible physics 
applications of these boundary conditions (or the corresponding delta function potential). We 
conclude in Sec. 6 with a brief discussion of the results and possible future directions. 

2 Edge states for the Dirac operator on flat manifolds with a 
hole 

In this section we give a couple of examples of the presence of edge states for the Dirac 
operator in the presence of boundaries along with analytical expression for their exact number 
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as function of the parameters of the problem. As briefly mentioned before, in the presence of 
boundaries, the self-adjointness issue for differential operators becomes non-trivial since one 
has to suitably define the domain of the operator (which includes specification of boundary 
conditions). For the Dirac operator the self-adjointness analysis in the presence of boundaries 
leads to the so called APS boundary condition mmm or its generalization thereof m- 
Since our aim is to illustrate the existence of edge states and in finding their exact number 
(as a function of the parameters in the problem) we consider situations where the system 
(including the boundary) has the isometries of a sphere. We first consider the problem on 
M 2 — V 2 , that is, on a two dimensional Euclidean spatial manifold with a disc removed. 
Next we consider the case of three dimensional spatial manifold with a ball removed (M 3 — 
B 3 ). Interestingly, we will find that whereas there is always at least one edge state in 2 + 1 
dimensions, this is not the case in 3 + 1 dimensions. 


2.1 Dirac operator on M 2 — V 2 

We are interested in finding the edge states for the Dirac operator on 2 + 1 dimensional 

Minkowski spacetime (with metric signature -I-) in the presence of a boundary on the 

spatial manifold in the form of a hole of radius rj, (that is, the spatial manifold is M 2 — T> 2 ). 
A suitable boundary condition will turn out to be the APS boundary condition mmm- 
The presence of edge states for Dirac operator on a disk with APS boundary condition has 
already been demonstrated in m and we will essentially follow the same analysis. In |13| 
only the existence of edge states was shown and more recently, in [28], a bound on these states 
was found. Below we give an exact counting of the number of edge states as a function of 
parameters appearing in the problem. 

The Dirac Hamiltonian in one particle quantum mechanics is 


H = — ia ■ V + mo/3, (1) 

where a 1 = 7 ° 7 * and /3 = 7 0 satisfy the condition (a *) 2 = /3 2 = 1 , i = (1,2) and we choose 
7 0 = <73, 7 * = ioi where ex’s are the Pauli matrices (note that our convention for the 7 ’s is 
different from that in m)- In polar coordinates the Hamiltonian is 

H = —ia r d r - — d e + m 0 a 3 , (2) 

r 

where a r = —02 cos 6 + 07 sin 6 and oq = 02 sin 6 + o\ cos 0 are explicitly given by 


( 0 ie~ ie \ ( 0 e~ ie \ 

Ur ~ V ~ie ie 0 ) ’ ae ~ \ e ie 0 ) ’ 


(3) 


It is easy to check that the total angular momentum operator J = L + S = —idg + ^a 3 
commutes with the Hamiltonian: [H, J] = 0, implying that the two can be simultaneously 
diagonalized. The (total) angular momentum eigenstates are: 


^ ~\e ) = e i{j -^ 6 



<//. + \e ) = 



(4) 


Expressing the eigenvectors 'Fj('r, 9) of the Hamiltonian ([ 2 ]) in terms of the angular mo¬ 
mentum eigenvectors Q: 


Vj(r,0) =Pj(r)^ + \e) + q j (r)i/)^ ] ( 6 ») 


(5) 
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leads to a pair of first order coupled equations for pj and qj: 

p'j ( r ) + ^ ^2 + mqj(r ) = E jqj (r), 

( 2 \ ^ f y , '\ 

-Qj(j') + j ~ x —"- mpj(r) = E jPj (r). 


(j + hf 


These can be decoupled to get the following two second order equations 

v 2 

Pj(r) = 0, 


Pj 


Pj ( r ) + y “ 


K - E)) + 


QUr) + 


Qj(r) 


( m o — + 






» = 0, 


( 6 ) 

(7) 

( 8 ) 

(9) 


where Ej is the eigenvalue corresponding to ^j(r,9). Equations ^ and © fall in the class 
of Bessel equations. We are interested in edge states or solutions that are localized at the 
boundary and these are obtained for e 2 = mg — Ej > 0 in which case the corresponding inde¬ 
pendent solutions are given in terms of the modified Bessel functions J-i (e 7 r) and K ■_ i (e,r). 
Demanding that the solution be square integrable on (?y,oo) leaves us with K-i(ejr) as the 
solution for edge states. 

An important issue when working on manifolds with boundaries is that of self-adjointness 
of the Hamiltonian operator and for this we will be following the APS prescription. With 
the bulk and the boundary inner products defined as (T,^) = and (ip,(p) = 

IdM dV^MV’V; respectively (here M = M 2 — V 2 and T |qm = ip and d *|dM = 4> on the 
boundary dM of M ), we need to evaluate — (AT,^) to find the condition for the 

self-adjointness of the Hamiltonian ([2]). 

Imposing single-valuedness of (\k, T) and using the fact that these functions fall-off to zero 
sufficiently fast as r —> oo we find that 

- (H®,®) =ir b J d0^(r b ,9)a r $(r b ,Q) = i{ip^(9),a r <p{9)), (10) 

where ip (9) = T(?y, 6) is the value of 'I'(r, 9) on the boundary at r = ry. For H to be 
self-adjoint, this should be equal to zero. 

Following the APS prescription, to make the Hamiltonian into a self-adjoint operator we 
need to find an operator AT(mo) on the boundary Hilbert space , H{dM) such that it anti¬ 
commutes with oy (which appears in the second term of the expression on right in (] 10[) ), 
K(m)a r = —oyA(m), such that the eigenvalue of A' 2 (mo) > 0. If such an operator exists 
then we can split T-L{dM.) = 'H\{dM)®'H 2 {dM) where T~L\{dM) and 'H. 2 {dM ) are orthogonal 
sub-spaces of 'H{dA4). 

Because K and oy anticommute, this would imply that for ipi{9) £ Wi(dA4), a r ipi(9) E 
'H 2 (dM). Thus, if we choose the boundary condition so that \k(7y, 9) = ip(9), 3>(ry, 9) = <p{9) 


are both in the same sub-space, say / H. 2 {dM), then the last expression in (10) will be equal 
to zero implying that the Hamiltonian is a symmetric operator (and also self-adjoint, as can 
be easily checked). 

It is easy to check that the operator 


K(m 0 ) = - 

n 


oy 


V ~ y ) + m 0 a 3 


( 11 ) 
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satisfies the first APS requirement in that it anti-commutes with o r . 

It has recently been suggested in m that a generalization of the APS boundary condition 
is possible and instead of working with the boundary operator K(mo) one can also work with 
a generalized boundary operator defined as 


K{n) = ia r K(n), 


( 12 ) 


where K{p) is as defined in © with parameter m o replaced by p. With this definition the 
explicit form of the K(p) operator is 


K(n) = —d e - 7T— + d^e- 
n 2 r b 


(13) 


It can be easily verified that this operator satisfies the first APS requirement {K(p), ay} = 0. 
However, note that K(p) and K(p) are related by a similarity transformation and therefore 
the main message of the generalization of m can be taken to be that one can replace the mass 
parameter mo in (11) by a more general parameter fi. In the following we will be working 
with the operator K(p) in (13). 

It is also easy to see that K(p) satisfies the other APS requirement K 2 (/i) > 0 as well, 
with the square of the corresponding eigenvalue given by 


\ 2 — + 


(14) 


The corresponding eigenvectors are 

e *O+i ) 0 


= Cj 


= d i 


i(j+h)o 




d \ 


rb 




+ 77 H^J'I ~ 7T 


rb 




, A j > 0 


, A j < 0, 


(15) 

(16) 


and the orthogonal sub-spaces H 1 (0A1) and can be identified as those corresponding 

to positive and negative eigenvalues respectively. 

Using the above expressions we now find the boundary conditions on pj and qj by equating 
the expression for \P / (r, 0 )\ r = rb as given in ([5]) to (16) (which corresponds to %~{dM) for 
K(p)): 


dj 




rb 


Pj(r b ) = dj, qj(r b ) = ( |Ay| + ( |Ay| + f- 


d 


n 


(17) 


This can be combined with the following relation between pj and qj (obtained directly from 
the Dirac equation) 


Pj{r) = - 


#’) , {j-Ddjir) 


+ 


(Ej + mo) (Ej + m 0 )r ’ 


to get 


N + ^ ) = 


p(Ej + mo) H- [j - - J { | Ay | + 


r b 


n 


Qj(,r b )- 


(18) 


The above boundary condition can be thought of as the generalization of the Robin bound¬ 
ary condition i^{r b ) = K^j'(r b ) with k in the present case being a function of the parameters mo, 
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Id and r b as well as a function of the spin j. After using the explicit form for qj (r) = cK ■ 1 (e^r) 

J 2 

(where c is a constant of integration) this can be written as 

fi(Ej + m 0 )K j _i(e j r b ) + ej (\Xj\ + K j+ i(ejr b ) = 0, (19) 

2 V r b J 2 

where, in writing the above expression we used the identity K' v (x) = vKAA — K v+ \{x). 


Counting the number of edge states 


We now count the number of edge states corresponding to the above boundary condition. We 


start by noting that since ( Ej + mo), ej and (|Ay| + j/r b ) are all greater than zero, (19) can 
be satisfied only for ju < 0. In other words, there will be edge state solutions only for fj, < 0. 


In (19) there appear three parameters (ii,mo,r b ) and it would seem that to get an exact 


analytic dependence of the number of edge states on all the three parameters might be difficult. 
However, by redefining (/r, mo) in terms of r b we can eliminate the radius r b of the hole from 
the problem. Specifically, we define 


M = 


r b 


m 0 = 


r b 


With this redefinition it is easy to verify that (19) can be written as 


i(±\Jb 2 - e j r b + b)K j _i(e j r b ) + ejr b (y/a 2 + j 2 + j)K j+ i(ejr b ) = 0, 


( 20 ) 


or writing ejr b = x and a = —A, A > 0 (recall that it was noted above that n < 0 to obtain 
edge states) we can rewrite the above equation as 


X (V A2 + j 2 + j)Kj + i(x) 

(±Vb 2 — x 2 + b)K-i(x) 
J 2 


( 21 ) 


In the above equation the ± sign in the denominator will correspond to j >0 and j < 0 
respectively. 

Next we note that the x-derivative of the rhs of the above equation in the limit x —> 0 
is zero (the only requirement for this argument to work is that j > 1/2 and j < —1/2 for 
the two cases mentioned above). Since the tangent to the rhs is zero as a function of x in 
the limit x —> 0 and that of lhs is obviously so (A being independent of x), the maximum 
number of bound state solutions will be given by finding the first curve (corresponding to the 
rhs of ©) which does not intersect the A = const, line in the limit x —> 0 (since all the 
curves corresponding to the rhs for different values of j start out horizontally by the above 
argument). 

We first consider the case j > 0 which corresponds to choosing the plus sign in the 
denominator of © (and as already mentioned, we need j > 1/2 for the argument to work). 
For x —> 0 the modified Bessel function can be approximated by 

lim K v {x) ~ f-') , y/v + \ > x > 0. (22) 

x —^0 2 \xJ 
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Using this on the rhs of (21) and remembering that j is a positive odd-half integer and that 


r(n) = (ri — 1)! we find that in the limit x —> 0 the condition in (21) implies 

Ab = (^j - 0 (j + y/A 2 + j 2 ). 


(23) 


The maximum number of allowed solutions j max will be given in terms of the solution of this 
equation 

A + b + 6a/1 + 4 A 2 + 8 Ab 

Jmax = 2(A + 26) ' ( j 


The other solution of (23) gives j < 0 and hence does not apply. 



Solution for the boundary condition (21) in the presence of generalized boundary operator in 


two dimensions for A = 3,b = 9, j > 0. 

Next we consider the case j < 0 corresponding to choosing the minus sign in the de¬ 


nominator of (21) (and as mentioned earlier, we need j < —1/2 for the argument to work). 


Since j < —1/2 we cannot use the approximation in (22). However, the modified Bessel 

functions satisfy the relation K v (x) = K- V (x) and therefore K-,i(x) = K_-_i(x) and 

J' 2 _ 2 


K-i(x) = K_,i (x). With these forms the condition in the approximation in (22) applies 
and we can now proceed to evaluate the rhs of the above equation in the limit x —>• 0. Here 
we also need to be careful about the factor in the denominator: b — (b 2 — x 2 ) 1 / 2 ~ x 2 /2b. 
With all these we finally get the equation 


- (j + \) A = Kj + Va*+P). 


(25) 


The solution of this equation which gives the minimum value of j (remember j < 0) that still 
leads to an edge state is 


Jmin — 


A + b + bV 1 + 4 A 2 + 8 Ab 
2(A + 2b) ’ 


( 26 ) 


which just differs by an overall sign from the corresponding expression in (24) for j > 0. 
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One can easily convince oneself that for j > 0 and when j max is not an integer, the actual 
number of edge states is given by the lowest integer greater than j max in (24). Similarly for 
j < 0 and when j m ; n is not an integer, the actual number of edge states is the greatest integer 
smaller than j m \ n in (26). When j max and j m j n are integers than the number of edge states 
is jmax — 1 and j m in + 1 since only for these values of j does the rhs of ( 21 ) intersects the 
A = const curve (the curves corresponding to j max and j m ; n are coincident with the A = const 
curve at x = 0 implying that there is no intersection). 

The correctness of the above counting (including the arguments leading upto it) can also 
be seen in Figs, [l] and [2] corresponding to j >0 and j < 0, respectively, for the parameter 
choice A = 3 and 6 = 9. From the figures we see that there are four edge states for each case. 
Our formulas ( |24| and p 6 | ) give j max = —Jmin = 3.69 and therefore our counting prescription 
implies that the number of edge states is 4 + 4 in agreement with what is seen in the figures. 



Figure 2: 


Solution for the boundary condition (21) in the presence of generalized boundary operator in 
two dimensions for A = 3,b = 9, j < 0. 


Number of edge states in various limits in 2 + 1 dimensions 


It might be of some interest to see the behavior of the number of edge states when the 
parameter A in (24) takes on certain limiting values (we only focus on j max , the case of j m ; n 
being identical). It is easy to see that 


• In the limit A —> 0, j max —> that is, their is only a single edge state. 

• In the limit A —> oo, j max —> b + b, that is, the number of edge states grows linearly 
with the mass of the Dirac particle when A is very large. 


2.2 Dirac operator on M 3 — B 3 

We next consider the Dirac operator in 3 + 1 dimensional Minkowski spacetime (with metric 

signature 4 -) with a three ball of radius r ^ removed from the the spatial manifold, that 

is, we work on the spatial slice M 3 — B 3 . Since the analysis of the Dirac equation in presence 
of spherical symmetry is standard, in the following we will be brief. The Dirac Hamiltonian 
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IS 


H = — ia ■ V + mo/3, 


with a 1 = 7 ° 7 * and /3 = 7°. We work in the representation 


— „,o 


a = 


0 (Ti 
Oi 0 


P = 


1 0 


0 -1 

Using the spherical polar coordinates the Hamiltonian can be written as 


ij 


*7* 


H = -if - -L-de - -!-d v + m. o7 u . 


(27) 

(28) 

(29) 

(30) 

0 -iip 

(31) 

Because of spherical symmetry the Hamiltonian commutes with the square of the total 
angular momentum J 2 , with = L + ( L being the orbital angular momentum and 

the spin), as well as with its ^-component J z , implying that these three operators can 
be diagonalized simultaneously. Following 
written as 


In the above equation 

r 

7 = 

with 


0 a r 
(J r 0 


1 ° = 


0 <70 

(70 0 


Y = 


0 < 7 V 
CTp 0 


(J 7* - 


cos 9 sin0e 1v 
sin 9e z,p — cos 9 


<70 = 


— sin 9 cos 9e 
cos 9e llf sin 9 


<7^0 — 


_J_ p up 
sin 0 c 


sin 0 


^1 (r,9,<p) = 


F( r ) y]™. 


the eigenvectors for the Hamiltonian can be 

GWY 


^2 (r,9,ip) =| 

-ig{r)y._ 1 

j 2 


(32) 


where 


y j r = y t 


i±|,m 


2 _ ( ± \J l±m + h Y ™ 5 


Y21 + 1 


lJfm+ \Y } 


(One can check that y J _i and y J , 1 are eigenstates of J 2 .) 

3 2 2 

Since we are interested in studying the effects of boundaries and boundary conditions, 
in the following we will concentrate only on d / i(r, 9, <p), the analysis for ^ 2 (r, 0 , 99 ) being 
analogous.The radial equations for F(r) and /(r) resulting from H^(r,9,f) = F\1 f (r,9,ip) 
are: 


(F - m 0 )F - 
(F + m 0 )/ + 


raf 

3> 

fif) ] 

_d/ + 

Ft 

r 

dF 


F(r)" 

dr 

V V 

r 


= 0 , 
= 0 . 


(33) 

(34) 


Eliminating f(r) in favor of F(r) and defining n = j — 1/2, m^ — E 2 = e 2 and e^r = x we get 
the equation 


( e i r ) 


,d 2 F(r) 

d(gjr ) 2 


+ 2e,r 


dF(r) 

d(ejr) 


— [{ejr) 2 + n(n + l)]F(r) = 0. 


(35) 
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This equation has two solutions c\I n+ i{ejr)/r 1 / 2 and C 2 -f£T n+ i (ejr)/r 1 ^ 2 of which only the 
second solution is finite on rj < r < oo (here c\ and C 2 are arbitrary constants). 

Following an analysis similar to that of the previous sub-section we find that the boundary 
operator which makes the Hamiltonian a self-adjoint operator is 


^o) = ^-^-^ + m o 7 0 , (36) 

n n r b 

or its generalization in the form K(fi) = i'y r K(p) (which is again related to I\{mo) by a 
similarity transformation) 


K(fi) 


7 r, , 7^ O , . -0 

— oq H-ci, + iii'f 

n r b 


1 

5 

n 


(37) 


where 7 ® = Yl 9 etc. We choose to work with the operator in 
condition resulting from it. 

The eigenvectors of K(/i) are 


(37) 


and find the boundary 


M + (M) 


M ( 0 , 4>) 



with the corresponding eigenvalues 





(38) 

(39) 


(40) 


As for the 2 + 1 dimensional case, if we now demand that '!'( r b , 9, tp) = ij)- (0, p) we end 
up with the boundary condition 


I'M H— 
r b 


J + l 


F\r b ) = 


n(E + m 0 ) + — (j-\ 
r b \ 2 


IM d- 

r b 


J + \ 


F(n). 


(41) 

From (35) we know that F(r) = cKj(ejr)/r l F. To simplify the subsequent expression(s) we 
now express p and mo in terms of r b as fj, = a/r b and mo = b/r b and define x = €jr b . With 
all this the equation for the boundary condition becomes 


—x 


a = 


\ja 2 + (j + 1) + (j + 1) 


. K j+ i(*) 


(Vb 2 — x 2 + b) 


Kj(x) 


(42) 


As in 2 + 1 dimensions it turns out that we get non-trivial solutions only for a < 0. Making 
the replacement a —> —A and comparing the resulting equation with the corresponding equa¬ 


tion ( 21 ) for 2 + 1 dimensions we find that the present equation can be obtained by replacing 


j with j + 1/2 in (21). All the previous arguments therefore go through and the number of 

+ 1/2 which gives the 


edge solutions can also be obtained by replacing j 
maximum number of edge states as 


in (24) with j D 


_-b + bV 1 + 4A 2 + 8 Ab 
Jmax “ 2 (A + 2b) ’ 


(43) 
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As an example consider A = 3, b = 9 for which the above equation gives j max = 3.1941. 
Graphically the situation is shown in figure [3j If we write j max = n m + /, where n m E Z and 
/ is the fractional part, then for / < 1/2, the number of edge states is n m and for / > 1/2 
the number of edge states is n m + 1. 



Solution for the boundary condition (42) in the presence of generalized boundary operator in 


three spatial dimensions for A = 3, b = 9. 


Number of edge states in various limits in 3 + 1 dimensions 


As for 2 + 1 dimensions, we can calculate the behavior of the number of edge states in various 
limits of the parameter A in (43) for 3 + 1 dimensions. We find that 


• In the limit A —> 0, j max —> 0. This is to be contrasted with the 2 + 1 dimensional case 
where one edge state survives in this limit. 

• In the limit A — > oo, j max —> b. This is similar to 2 + 1 dimensions where the behavior 
in this limit is b + 1/2. 


3 Singular potentials and the boundary effects 

In this section we suggest an interesting possibility whereby the effect of the boundary and 
the corresponding boundary condition can be mimicked by replacing the boundary and the 
boundary condition by delta function potential with an appropriately chosen strength at the 
location of the original boundary. In essence, what is being suggested is that if the original 
domain for the operator was restricted to the interval n, < r < oo with suitable boundary 
condition at r = then one can extend the domain to the interval 0 < r < oo (if we 
think of r as the radial coordinate) and introduce a delta function potential P5(r — r/,) in the 
Hamiltonian with a suitably chosen value for P. This has the additional advantage that in the 
absence of a boundary it becomes easier to address the self-adjointness issue. From physics 
perspective also this makes sense since, in general, boundaries (and boundary conditions) are 
achieved by introducing suitable potentials in the experimental set up. 


12 








To set things in context we first recall the text book example involving the Laplacian in 
one dimension where, in presence of the delta function potential one can get a single bound 
state. We will show that in higher dimensions the number of bound states can depend on 
the strength of the delta potential. As has been shown in [13] it is not necessary to restrict 
to only 5 potential and one can also include its derivative. The key point is that instead of 
working with the Schroedinger equation, one uses the conservation of the probability current 
(in other words, one uses the continuity equation). Interestingly, while the <5 function allows 
for the discontinuity in the derivative of the wave function the 5' potential allows for the 
discontinuity in wave function itself. But the probability remains continuous. Although in 
the context of the Laplacian the inclusion of the 5' term might appear a bit adhoc, we will 
see that for the Dirac operator the 5' term arises naturally. 

To have a concrete quantitative understanding of what is being proposed we first consider 
the case of the Laplacian on M 2 — D 2 . This problem was considered in [9j where an exact 
expression was obtained for the number of bound (edge) states in the presence of the Robin 
boundary condition Kip+d^tj} = 0. We now consider the problem where the disc T> 2 is restored 
and the Robin boundary condition is replaced by a delta function potential V(r) = P5(r — rb ) 
where P is the strength of the potential. 

The Schroedinger equation takes the form 


d 2 i/j 1 d 2 ip 

Q r 2 r 2 QQ2 

For r / r b this equation is solved by 


1 dib 

H- 7 T~ + PS(r — 7 , b)'i/> 

r or 


A 4). 


ip(r<r b ,6) = cl n (\/\r)e m9 , 
i/j(r>rb,9) = dK n (\/Xr)e ine . 


(44) 


(45) 

(46) 


At r = r b we need to worry about the delta function singularity and we integrate the radial 
equation (obtained after using if) = R(r)e in9 in (44)) from r b — e to r b + e and take the limit 
e —> 0. This gives the equation 


d R 
dr 


r b +£ 


+ PR(r b ) = 0. 


(47) 


J r b -e 


In writing this equation we have assumed the continuity of the wave function across r = r b 
which implies cl n {\f\rb) = dK n (\f\rb). Using I' v {x) = I v -i{x) — vl v {x)/x and K' u (x) = 
—K u _i(x) — vK v {x)/x (where prime implies derivative with respect to the argument) then 
gives the condition 


Pn— 1 (^) . In— 1 (■£*) 
K n {x) I n (x) 


(48) 


where we have defined P = Q/r b and x = \f\rb- 

Now it turns out that the tangent to the curves corresponding to the expression on the 
left of the above equation is horizontal near x = 0 and therefore the arguments used in 
the previous section to find analytic expression for the number of edgestates can be applied. 
Using the approximation in (22) (and a similar approximation for I u (x)) we find that the 
above equation in the limit x —> 0 leads to the following expression for the maximum number 
of bound states: 


lim x 


x 


2 (rim - 1) 


+ 


2 n n 


— 2'ZT'm — Q • 


(49) 
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On the other hand, in [9], in presence of the Robin boundary condition the maximum 
number of edge states was found to be n m = — (k < 0). Comparing this with the 
previous expression we find that the problem in the presence of delta potential but with no 
boundaries is equivalent to the problem in the presence of Robin boundary condition if we 
choose Q = —2K,rb or equivalently 

P = -2k. (50) 

Having established the correspondence between Robin boundary condition and the pres¬ 
ence of a delta function potential (of appropriate strength) at the location of the boundary 
for the Laplacian, the next task would be to do the same for the Dirac operator with APS 
boundary condition. However, before doing that we would like to point to the analysis by 
Boya and Sudarshan M where it is shown that the usual matching conditions employed in 
problems involving delta potential can be seen as arising naturally from the current conser¬ 
vation equation (the continuity equation) for the Schroedinger equation. Additionally, they 
show that use of current conservation as the basic principle allows for more general singular 
potentials - those involving both 5(r) and 8'(r) type singularity - and these can be handled 
by essentially similar methods. Inclusion of 5'(r) type potential directly at the level of the 
Schroedinger equation, requires careful analysis. 

From the point of view of the Laplacian or the Schroedinger equation, inclusion of a 5'(r) 
potential might seem a bit adhoc. However, we next show that when considering the Dirac 
operator the S'(r) type term in the potential arises naturally. We begin by considering the 
first example of the previous section but now with the hole of radius r {, replaced by a delta 
function potential V(r) = —P5(r — r&). The Hamiltonian is therefore the same as in © 
except for the addition of V(r): 

H = —ia r d r - -d$ + mas + V(r). (51) 

r 

As usual, we want to solve the equation ) = Ejipj(r,0) and since the system is 

still spherically symmetric, the wavefunction is still expressed as in (pi). In the presence of a 
potential V(r) = P5(r — r^) the first order Dirac equations ([b]) and ([7]) are modified by the 
addition of the potential term to 

Pj (r) + (j + 0 + mqj (r) + V(r)qj (r) = E jqj ( r ), (52) 

~q'j(r) + (j - 0 - mpj (r) + V(r) Pj (r) = E jPj {r). (53) 


These equations can be decoupled as before and, for instance, leads to the following second 
order equation for qj (r) (compare with ([9]) corresponding to V(r) = 0) 


q"(r) + 


q'j( r ) 


(m 2 - (Ej - P5(r - n)) 2 ) + 


rq'j(r) - ( j - - qj{r) 


P5'{r — r&) 
r(Ej — V(r) + m) 


C7 - h) 


= 0 . 


D 2 n 


qj( j 


(54) 


Since the operator above is essentially the Laplacian, we see that inclusion of the delta poten¬ 
tial in the Dirac operator automatically leads to the presence of 5'(r) terms for the Laplacian. 
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We can now go on to show that even for the Dirac operator the effect of the APS boundary 
condition, as discussed in the previous section, can be mimicked by the delta potential. 
In other words, as for the Laplacian, we can choose the strength P of the delta potential 
appropriately so that the number of edge states for the two problems are identical. For r < r b 
and for r > r b , V(r) = 0 and, therefore, in these regions the second order equations for pj 
and qj will be the same as in ([8]) and ([9]) with the solutions given in terms of the modified 
Bessel functions. Specifically, if we choose the (Jj(r') solution as 

Qj(r < r b ) = cil-i(ejr), qj(r > r b ) = c 2 K-i(ejr), (55) 

J 2 J 2 

then the corresponding solutions for pj (r ) (using 0 ) are 

Pj(r < r b ) = - Cl , Cj = Jj + i(e 3 -r), Pj (r > r b ) = - C * 3 K j+ i(e jr ), (56) 

mo ± J rriQ — ej m 0 ± Jmfi-ej 


ml - E 2 -. 


where €j = 

As for the Laplacian (Schroedinger equation), we might try and integrate (54). However, 
the presence of the 5' terms does not allow a straight forward use of this standard procedure. 
One might, presumably, generalize the method of Boya and Sudarshan [14] and use the 
current conservation equation for the Dirac Hamiltonian. However, in this paper we follow 
the procedure of [30l ED and directly work with the first order equations ( |52[ ) and ( |53[ ). 

We multiply (52) by pj and (53) by qj and add the two equations to obtain: 


W , (ffl «]) (p? + # - 

2 + 2 r 2 r PjQj ' 

We integrate the above equation from — e to + e and take the limit e 

2,2 2,2 

Pj> + Qj> = Pj< + Qj<, 


(57) 

0 to the condition 

(58) 


where pj ^ = Pj (rb e) and qj ^ = qj (rb 4= e) ■ This condition simply fixes the normalization of 
the solution on one side r = in terms of the normalization chosen on the other side. 

Next we multiply (52) by qj and (53) by pj and subtract the two equations. This gives: 

PjQj ~ PjQj = 2// ; V/? - (Ej ~ m)q 2 - (Ej + m)p 2 - P6(r - r b )(p 2 + q 2 ). (59) 

We divide the above equation by p 2 + q 2 and, as before, integrate the resulting expression 
from r b — e to r& + e. On the rhs, except for the last term all other terms give zero, and we 
get the second condition 


tan 


-l 


l n+t 


Op 

-Pj J r b -e 


= ~P , 


(60) 


or defining P = tan 1 Q, the above equation can be rewritten as 


Qj> 

Pj> 


Qj< 

Pj< 


— Q ) (i + Q 


> Qj< 
'Pj< 


-i 


(61) 
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Using the explicit solutions (55) and (56) and (61) in the e 
can be rewritten as 


0 limit, the above condition 


Q = - 


(b ± \Jb 2 — x 2 ) 






+ 


a 


i+sw 


i - 


(i6 ± \/6 2 - .t 2 ) 2 (*) ( x )' 

z 2 / i+ i(s)A i+ i(®) 


-l 


(62 ) 

Note that in this equation the normalization factors ci and C 2 do not appear anywhere as only 
the ratio of the solutions are involved. As in the previous section, we have used the definition 
?no = b/rb and eyri, = x. 

Next we need to find the edge states from the above equation. As before, we find that 
the slope of the curves corresponding to the rhs of (62) in the limit x —> 0 is zero and, 


therefore, the same procedure can be used to find the number of edge states (and, as before, 
we only require / > 1/2 for the argument to work). We start by noting that the two signs 
in — x 2 correspond to positive and negative energy solutions. We therefore have four 

possibilities - (i) Ej > 0, j > 0 (ii) Ej > 0, / <0 (iii) Ej < 0, / >0 and (iv) Ej < 0, / < 0. 
Following the earlier method we find that these four possibilities will lead to two solutions 
each for j max and j m ; n . For instance, the first case above where we choose the positive square 
root and j > 0, the total number of bound states is found to be 


Jmax,l 


— Qb + 


1 


(63) 


Similarly, case two leads to 


Jmin,l — Qb . 


(64) 


Consistency with the requirement that jmax ,l > 0 and j mm ,i < 0 implies that Q > 0. In a 
similar manner the other two cases, corresponding to choosing the minus sign for the square 
root, lead to j max ,2 = —Qb — 1/2 and j m i n ,2 = Qb— 1/2 and consistency implies that for these 
two cases Q < 0. 

Now the sign of Q (equivalently of P) determines whether the particle sees a potential 
hill or a potential well. However, this interpretation also depends on whether the particle 
has positive energy or negative energy. The first two cases corresponding to positive energy 
and Q > 0, therefore, imply that the particle sees a potential well. Similarly, the other two 
solutions for Q < 0 correspond to a negative energy particle seeing a potential well. Thus we 

find that /max,l = /min,2 and /max,2 = /min,l- 

If we now equate the expression for one of the / maX ) say / max , i in (63), to the number 


found in (24) in the presence of the APS boundary condition, we get the condition 


Vl + 4A 2 + 8 Ab 
Q ~ 2 (A + 2b) ’ 


(65) 


which leads to the same number of edge states with / > 0 for the problem with the delta 
potential as obtained in presence of the APS boundary condition. And as shown above, the 
negative of this Q value leads to the same / m in ,2 as found with the APS boundary condition. 
It is interesting to see that the role of the sign of the potential is seen only for the edge 
states. The expression for Q will be slightly modified for /max ,2 (and / m in,i) and is given by 
Q = {2 A + 3b + bQl + 4 A 2 + 8 Ab)/2b(A + 26). 
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4 Time dependent (moving) boundaries and boundary condi¬ 
tions 

The problem of quantum mechanics and quantum field theory on manifolds with time chang¬ 
ing or moving boundaries is interesting in its own right but also has important physics impli¬ 
cations. For instance, it might be of interest to study Casimir effect in the presence of moving 
boundaries and, in the context of general relativity and cosmology also such situations arise 
naturally in the form of evolving horizons (black hole or cosmological) which are often treated 
as boundaries. Moving mirror and associated Unruh effect are other examples of phenomena 
involving moving boundaries. In this section we, therefore, consider the question as to what 
happens to quantum mechanical problems when the boundary is changing with time? This 
change can be in variety of ways. There can be expansion, translation or there can be a 
change of shape of the boundary. 

The difficulty with quantum mechanical problems involving moving boundary is that the 
domain of the Hamiltonian operator in the space of £ 2 (M) (in order that the Hamiltonian 
be a self-adjoint operator) changes with time. As the boundary evolves, the wave function 
will go outside the domain and hence the description of time evolution becomes subtle. To 
deal with this we have to deal with the whole Hilbert space and the expansion of functions 
in terms eigenfunctions will change with time (26j . 

4.1 Expanding disc 

Below we consider the situation involving M 2 — T> 2 with the radius of the disc being a function 
of time and explore how to do quantum physics in such a situation. We focus on the case 
where the Hamiltonian is Laplacian in the region M 2 — D 2 and zero outside this domain (from 
the analysis of the previous sections it should be clear that situation involving the Dirac 
operator is not fundamentally different from that of the Laplacian). 

The expanding radius is given by 

R(t) = e?®Ro, /(0) = 0, ( 66 ) 

and we will be considering imposition of the Robin boundary condition with parameter k. 
The basic idea would be to translate the problem to one involving fixed boundary but a time 
dependent Hamiltonian. For this consider the following transformation of the wave function 

X(r,0) = Uip{r,9), (67) 

where U(r,9 ) = e~ a ^ r \ Then 

x\dD = e~ a i>(r, 9 )\q D . (68) 

On the other hand, the normal derivative of x is (remember that we are imposing Robin 
boundary condition) 

dnX\dD = -d n a\ dD (e~ a ^(r, 9)) \ dD - e~ a(r) 9)\ 9D = -{k + d n a)x\dD , (69) 

where we have made use of ijj'(r)\QD = — Kip(r)\gD- Thus we see that the transformation 
effects a change in the boundary condition (and thereby the domain in the Hilbert space) by 
changing k —> k 1 = n + d n a\gr) (note that the prime on k does not denote derivative with 
respect to r). 
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Now the above transformation implies that the Hamiltonian changes to H' = UHU _1 . 
Using this we can maintain the domain of the Hilbert space same while changing the Hamil¬ 
tonian and hence the spectrum. If the a depends on time the Hamiltonian will change to 

r)TJ 

H’ = UHU~ 1 + i—U~ 1 . (70) 

It will be a time dependent Hamiltonian. 

We want that k' , the new parameter for the Robin boundary condition when the problem 
is translated to one with a fixed boundary i?(0), be such that the physics for the problem 
(«/, -R(O)) is the same as that for (k, R(t)). In particular, this means that the number of edge 
states for the two problems should be the same. As was shown in [9], the number of edge 
states for the Laplacian on M 2 — V 2 is — nr^, {jr^ being the radius of the boundary) and this 
implies that 

k'R{ 0) = nR(t). (71) 

Hence n' = = e? k. But we also know: 


d n a = k — k. 


(72) 


Hence we get 

a = (k' — n)r = (e? — 1 )nr. (73) 

Having obtained a(r,t ) we can write U explicitly and hence the time dependent H'. The 
transition amplitudes are the important quantities that can be found using the Hamiltonian 
which is time dependent. 


4.2 Delta function potential 

It is well known that in one dimension the boundary conditions can be obtained from a delta 
function potential and, in the previous section, we already demonstrated that the same holds 
even in higher dimensions. In principle the strength of the potential can vary along the 
boundary circle so that boundary conditions have angular dependence. We will ignore this 
possibility and take the strength to be uniform along the boundary circle. 

Consider the following Hamiltonian in M 2 

H = -A + g5(r- R{t)). (74) 

Here the delta function along a curve C is defined as: 

J 5(f(C))'i/j(x)d 2 x = j ip(x)dx. (75) 

For example when i/j(x) = 1 and curve C is a circle of radius R we get 2nR. Now R(t) = 
R( 0) and if we apply dilation operator D such that r —> kr then we get for the potential 
6(hr — kR(t )). If we fix k to be the potential becomes 5(e~^r — R(0)) = e^5(r — R( 0)). 
The time dependence is shifted to the strength of potential. This is analogous to changing 
the Hamiltonian to a time dependent one by keeping the domain of the Hilbert space same 
for all times. 
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4.3 Solvable time dependent boundary problems 

Interestingly, the case of one dimensional delta function potential whose strength varies with 
time can be exactly solved in several cases. Several of these time dependent potentials cor¬ 
respond to problems with moving boundaries, for example, where the boundary moves with 
uniform velocity or has uniform acceleration. In the first case the Galilean invariance pro¬ 
vides the solution. In the second case the Hamiltonian can be transformed to a constant 
gravitational field whose solution is given by Airy functions [261123|. 

There are several examples in higher dimensions where the time dependent problem can 
be exactly solved which corresponds to moving boundaries. One such example is given by 
Berry and Klein [23] : 

H (r,pJ(t)) = ^ n + ^ v (r/ 1 ). (76) 

Here l depends on time. This Hamiltonian is not conserved. If the time dependence is of the 
form /(f) = V at 2 + 2 bt + c then one can write in a comoving frame 

H(p,n,k) = +V{p) + ^kp 2 , (77) 

where p = r/l and k = m(ac — b 2 ) which is conserved in p,r = The expanding 

disc in M 2 and ball in M 3 will come under this class of Hamiltonians. A detailed analysis 
of the quantum mechanics with implications for field theory as well as extension to Dirac 
Hamiltonian will be presented elsewhere. 

5 Applications 

Manifolds with boundary is considered in several applications. We will list few of them. 
Analysis in several new situations will be presented elsewhere. 

• Quantum Hall effect: This is a classic example of electron gas in a two dimensional 
manifold with boundary. Naturally the edge excitations are chiral and has been analysed 
to a large degree. Topological insulators arising in 2d and 3d also can be understood 
in this context. See for example m- But moving boundaries in this context have not 
been considered both theoretically and experimentally. 

• Casimir effect: Casimir effect of force between plates due to vacuum fluctuations can be 
analysed using self-adjoint extensions of the Laplacian or Dirac operators with suitable 
boundary conditions. Effects of modelling the plates in terms of singular potentials 
has been investigated in several recent papers Hang nans Eg. Asorey et al [321 
have analyzed the effect of boundary conditions on the Casimir force and have showed 
that depending on the nature of the boundary condition this force can be attractive, 
repulsive or it can also vanish. The focus is on positive Laplacian operator and it would 
be interesting to see the effect of finite temperature and self-interactions. Detailed 
changes in the nature of the boundary like half-line or boundary with holes have been 
considered in the literature [33] . The nature of this effect with moving boundaries will 
be new tool for testing the developments posed here. 

• Isolated and cosmological horizons: Blackholes present novel and unstable behaviour 
classically due to the appearance of null surfaces such as horizons. They require quantum 
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theory to understand several questions. Membrane paradigm |34[l35j is one framework in 
which the horizon acts as boundary and requires boundary conditions. One can analyse 
the situation using self-adjoint extensions for suitable boundary conditions. Similar 
situation arises for the cosmological horizon of de Sitter space of evolving universe. The 
nature of evolving boundaries in this situation can be used to relate the stationary 
coordinates and ingoing observers. 

• Edge excitations localised in the boundary have dramatic effects in gas of bosons as 
pointed out by one of us [36] • Evolving boundaries in this case and effects on statistical 
mechanics of gas of bosons/fermions require further analysis. A nice example of edge 
states actually appearing and requiring clear understanding of boundary conditions is 
provided by plasmons [371 ESI [39]. 


6 Discussion 


Problems requiring solution of quantum mechanical equations on manifolds with boundaries 
feature very often in physics. The most important issue in these situations is that of the self¬ 
adjointness of the corresponding Hamiltonian operator on the given domain. Although the 
two most commonly employed boundary conditions in the physics literature are the Dirichlet 
and the Neumann boundary conditions, it is known that self-adjointness of the Hamiltonian 
permits more general boundary conditions. For the case of the Laplacian these are the Robin 
boundary conditions while for the Dirac operator these are the APS boundary conditions. 

In this paper we showed that the APS boundary condition applicable for the Dirac operator 
can be thought of as generalized Robin boundary condition. The Robin boundary condition 
is usually given in the form ip'(r)\oD = — ni^{r)\dD where k is a dimensionful parameter and 


is usually taken to be a constant. If we now look back at the equation (18) which implements 


the APS boundary condition in 2 + 1 dimensions (and the corresponding equation (41) for 
3 + 1 dimensions) then it is obvious that these are essentially of the same form as the Robin 
boundary condition. 

However, there is an essential difference. While in the usual application of the Robin 
boundary condition the parameter k is taken to be a constant, we saw that, more generally, the 
parameter can depend on the spin (angular momentum) state j as well on the corresponding 
energy eigenvalue in the given state (in addition to the mass parameters mo and /r). We gave 
an exact counting for the number of edge states for the Dirac operator in presence of the APS 
boundary condition in 2 + 1 dimensions as well as in 3 + 1 dimensions for the case when the 
boundary has the isometry of a sphere. 

It is known that the issue of self-adjointness of operators becomes much more intricate in 
the presence of boundaries. One also expects that in actual physical situations the boundaries 
and boundary conditions would be implemented by suitable interactions. Based on this idea 
we showed that the effect of boundary conditions can be mimicked by extending the spatial 
manifold beyond the boundary (so that there is no boundary any more) and by putting a 
delta function potential of suitable strength at the location of the boundary. We worked out 
this correspondence explicitly for both, the Laplacian and the Dirac operator. For the case of 
the Laplacian we showed how the strength of the delta potential is related to the parameter 
k appearing in the Robin boundary condition, so that both the problems lead to the same 
number of edge states. 


20 




For the case of the Dirac operator the correspondence is essentially there except for one 
caveat that while in the case of the APS boundary condition the number of edge states with 
j > 0 is the same as the number of states with j < 0 (for a given choice of the parameters), 
for the case where the boundary condition is replaced by a delta potential, the number of 
edge states (for a given value of the strength of the delta potential) for j > 0 is one more 
than that for j < 0. Here it is also useful to note that like the parameter k of the Robin 
boundary condition, the strength of the delta potential is a dimensionful quantity in the case 
of the Laplacian (a second order operator) whereas for the Dirac operator, the strength of 
the delta potential is dimensionless since the Dirac operator is a first order operator. 

Finally we showed how the complicated problem of moving boundaries can be translated 
to a different form. In the presence of moving boundaries, the domain of the Hamiltonian 
changes with time and therefore the question of self-adjointness of operators becomes much 
more complicated [26j. We showed that it is possible to make a transformation such that the 
boundary becomes fixed and the effect of the moving boundary is contained in a modified 
Hamiltonian and boundary condition. 

As previously alluded to, their are numerous situations where one is required to do quan¬ 
tum mechanical calculations on manifolds with boundaries (including moving boundaries). 
These include quantum Hall systems and topological insulators on the condensed matter side 
and cosmological and black hole horizon on the gravity side. For instance, use of self-adjoint 
extensions for the Dirac Hamiltonian on graphene with defects was considered in [30]. Addi¬ 
tionally, interesting question regarding heat kernel expansion for scalar field theory on finite 
interval has been analyzed in m following the earlier work in [42] . As in [2T1 22] , only posi¬ 
tive Laplacian is considered to obtain heat kernel coefficients. It will be interesting to extend 
this analysis to finite temperature condition or with self interactions. In a future work we 
hope to investigate the implications of the analysis of this paper on some of these physical 
systems. 
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